Unipolar transport and shot noise in metal-semiconductor-metal structures 
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We carry out a self-consistent analytical theory of unipolar current and noise properties of metal- 
semiconductor-metal structures made of highly resistive semiconductors in the presence of an applied 
bias of arbitrary strength. By including the effects of the diffusion current we succeed to study the 
whole range of carrier injection conditions going from low level injection, where the structure behaves 
as a linear resistor, to high level injection, where the structure behaves as a space charge limited 
diode. We show that these structures display shot noise at the highest voltages. Remarkably the 
crossover from Nyquist noise to shot noise exhibits a complicate behavior with increasing current 
where an initial square root dependence (double thermal noise) is followed by a cubic power law. 



I. INTRODUCTION 

Metal-semiconductor-metal structures have been of 
fundamental and applied interest since the birth of the 
physics of semiconductor device. Here we concentrate on 
the relevant case of high resistivity semiconductors where 
metals usually form Ohmic or injecting contacts.^ A com- 
mon feature of these structures is that carrier transport 
is mainly limited by diffusion, with thermionic emission 
playing a negligible role. Furthermore, when a bias volt- 
age is applied to these structures, they can behave either 
as a linear resistor (if carrier injection from the contacts 
is negligibly small) or as non-linear space charge limited 
device (if carrier injection from the contacts is extremely 
high).^ Therefore, these structures are still offering ex- 
cellent opportunities to study current noise properties 
induced by diffusion noise in a variety of transport con- 
ditions. 

The present knowledge of the noise properties of these 
structures can be summarized as follows. At thermody- 
namic equilibrium, the low frequency spectral density of 
current fluctuations, Sj''{0), is given by 



5f (0) - 



(1) 



in agreement with Nyquist theorem."^ Here, R'^'^ is the 
equilibrium device resistance, T the bath temperature 
and ks the Boltzmann constant. 

Out of equilibrium, when an external bias (voltage or 
current) is applied, Nyquist relation no longer holds and 
deviations from Eq.(l) are expected. These deviations 
have been studied in detail in two limiting cases, namely: 
(i) when the structure behaves as a space charge limited 
diode (with a strong inhomogeneous profile of the free 
carrier density), '^"^ and (ii) when the structure behaves 



as a linear resistor (with a homogeneous profile of the 
free carrier density).^ 

In the former case (i) the current / displays a quadratic 
dependence on the applied voltage V (Mott-Gurney law) 
/ — PV"^, where (3 is a sample dependent parameter. Ac- 
cordingly, the low frequency spectral density of current 
fluctuations was found to take the form^~^ 



5/(0) = 



Z(0) 



1/2 7-1/2 



(2) 



where Z{Q) is the low frequency impedance (differential 
resistance) . From the form of the first equality, this type 
of noise was called double thermal noise.* The cross-over 
between Nyquist noise and high voltage space charge lim- 
ited conditions is summarized by the formula 



5/(0) = 4fcBr 



V 



dl_ 

dV 



(3) 



Highly resistive structures under strong carrier injection 
conditions have been used to successfully test experimen- 
tally this prediction.* 

In the latter case (ii) the structure behaves as a linear 
resistor and displays a resistance R — L/(qAfin), where 
A is the cross sectional area, q the carrier charge, /i the 
mobility and n the free carrier density (the bar denotes 
an average with respect to fluctuations). When the ap- 
plied bias is high enough as to make the transit time due 
to drift tt = L?' / fiV significantly shorter than the dielec- 
tric relaxation time = e/qfj/n, where e is the dielectric 
permittivity, the structure was predicted to exhibit shot 



with 



5/(0) = 2ql. 



(4) 



The above results show two completely different be- 
haviors for the non-equilibrium current noise proper- 
ties of metal-semiconductor-metal structures depending 
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on whether the stationary free carrier distribution is 
strongly inhomogcncous (double thermal noise) or homo- 
geneous (shot noise). The physical reason for such a dis- 
crepancy of behavior is an intriguing problem at present. 
In the authors opinion its origin can be traced back to 
the theoretical difficulty of accounting for the effects of 
the diffusion current on the non-equilibrium current fluc- 
tuations in the whole range of different carrier injection 
conditions. 

The purpose of the present paper is to overcome 
such a difliculty by developing an analytical theory for 
the non-equilibrium current noise properties of mctal- 
semiconductor-metal structures made of highly resistive 
semiconductors that includes properly the eflfects of the 
diffusion current. The theory is valid in the whole range 
of physical conditions going from homogeneous to highly 
inhomogeneous free carrier distributions, thus allowing 
to describe in a continuous manner the current spectral 
density as a function of the applied bias for different lev- 
els of current injection. 

The paper is organized as follows. Section II describes 
the system under study. Section III presents the physi- 
cal model used to describe the low frequency transport 
and the non-equilibrium current fluctuation properties. 
Section IV is devoted to analyze the charge transport 
properties of the device, while Section V analyzes its non- 
equilibrium current noise. In particular, we analyze the 
effects of space charge on the shot noise properties. Sec- 
tion VI draws the main conclusions of our investigation. 
Two appendices provide the most technical derivations. 



II. SYSTEM UNDER STUDY 

The system under study consists of an active region 
made by a non-degenerate n-type semiconducting ma- 
terial sandwiched between two metal plates that act as 
contacts. The semiconductor is assiimed to be lightly 
doped with a donor density Nd- All donors are assumed 
to be ionized at the considered temperatures (the case of 
a p-type semiconductor can be analogously considered) . 
To justify a one-dimensional electrostatic treatment in 
the X direction and to neglect the effects of the bound- 
aries in the y and z directions the transversal size of the 
sample is assumed to be much larger than the character- 
istic Debye screening length. The metals are assumed to 
form ohmic injecting contacts so that the voltage drop in- 
side them, or equivalently, the contact resistance, can be 
neglected to a good approximation. Accordingly, when a 
voltage is applied to the structure, all the potential drop 
takes place inside the semiconductor between positions 
X = and x = L, and the contacts can be excluded 
from consideration. Depending on the parameter values 
of the contacts homogeneous as well as inhomogeneous 
conditions will be studied. 

We assume that inelastic scattering processes with 
phonons are dominant, so that carrier thermalization at 



the bath temperature holds at any point. Accordingly, 
for the applied bias considered here no carrier heating 
takes place which allows us to use a fleld independent 
electron mobility and diffusivity. 



III. PHYSICAL MODEL 

A. Charge transport 

The transport approach appropriate to describe the 
electrical properties of the metal-semiconductor-mctal 
structure under study consists of the standard drift- 
diffusion current equation self-consistently coupled to 
the Poisson equation and supplemented by appropriate 
boundary conditions. The current drift-diffusion equa- 
tion reads^ 



- = qn{x)nE{x) + QD—^, 



(5) 



where A is the cross sectional area, n{x) the local car- 
rier density, fjL the mobility (assumed to be field indepen- 
dent), -E'(.x) the local electric field and D the diffusion 
coefficient (related to the mobility through Einstein's re- 
lation D/ n = ksT/q, where non-degenerate statistics is 
assumed). The bar denotes an average with respect to 
fluctuations. The Poisson equation is as usual 



dE{x) 
dx 



■ n{x)] 



(6) 



with e the static dielectric constant of the semiconduct- 
ing material. The appropriate boundary conditions to 
describe the ohmic injecting contacts are given by (see 
appendix A) 



n(0) = n{L) = ric 



(7) 



where the carrier density at the contacts ric is indepen- 
dent of the applied bias and given by (see appendix A) 



Tin 



Nc exp - 



(8) 



Here, Nc is the effective density of states in the con- 
duction band and 0t,„ the metal-semiconductor barrier 
height. For values of the contact parameters such as 
ric = the stationary free carrier distribution is ho- 
mogeneous, while for values such as ric > Nd there is 
net carrier injection from the contacts and the station- 
ary carrier distribution is inhomogeneous. 

Equations (5) and (6) can be combined into a single 
equation for the electric field 



-D 



<PE{x) - 



dx^ 



- iiE{x) 



dE{x) qNo 
dx e 



subject to the boundary conditions 
dE{0) _ dE{L) 



dx 



dx 



= -A 



(10) 
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B. Current fluctuations 

The only source of fluctuations in the system is re- 
lated to the difFusion of carriers inside the structure. Ac- 
cordingly, the low frequency noise properties (beyond 
1// noise) can be described through a drift- diffusion- 
Langevin model. It consists of the linearized version 
of the transport model presented in the previous subsec- 
tion supplemented by the appropriate Langevin source 
which describes the diffusion noise. The current equa- 
tion in explicit form reads,®"^^ 



SI{t) 



qfiE{x)6nx{t) + qn{x)n8Ex{t) + 



qD 



dx 



A 



(11) 



where SE^ (t) and Srix {t) refer to the fluctuations of elec- 
tric field and number density at point x, respectively, and 
51{t) refers to the fluctuations of the total current. In 
Eq.(ll), consistently with the low- frequency limit taken 
here, the displacement current is neglected. Moreover, 
SIx{t) is a Langevin source associated with the fluctua- 
tions of the current induced by the diflfusion of carriers 
inside the sample. It has zero mean and low frequency 
spectral density, 



/ + 00 
-OO 



dt {5h{t)5I^> it')) = K{x)5{x - x'), (12) 



with brackets indicating ensemble average, and 
K{x) = AqAksTimix) = 4:Aq'^Dn{x), 



(13) 



where in the last equality use is made of the Einstein re- 
lation. Finally, the linearized Poisson equation is given 

by 



— = — bnx{t), 

dx e 



(14) 



and the linearized boundary conditions by (see appendix 
A) 



driLit) = Snoit) = 0. 



(15) 



In analogy with the case of the transport equations, 
we can combine Eqs.(ll) and (14) into a single equation 
for the electric fleld fluctuation of the form 



SI{t) - 6Ut) ^ ^ dHE^jt) _ ^-^^^^dSE^it) 



dx^ 
-fxSE^it) 



dx 



dE{x) qNo 
dx e 



with boundary conditions 



dSE^it) 



dx 



dSE^t) 



dx 



0. 



(16) 



(17) 



IV. STATIONARY SPATIAL PROFILES AND 
CURRENT- VOLTAGE CHARACTERISTICS 

In this section we present the results concerning the 
transport properties for the structure under study calcu- 
lated from the model presented in Sec. Ill A. In particu- 
lar we will focus on the transition from homogeneous to 
inhomogeneous conditions. To this end, for given prop- 
erties of the semiconductor we will vary the contact den- 
sity from ric = Njj (corresponding to an homogeneous 
condition with Ohmic contacts) to values Uc > (cor- 
responding to inhomogeneous conditions with injecting 
contacts). Note, that for given properties of the semi- 
conductor, to vary ric corresponds to vary the contact 
barrier height </>(,„, i.e. to change the material for the 
metal contact. 

In what follows, the results are conveniently discussed 
in terms of the values of two dimensionless parameters 



/ = L/Ln and a = Uc/Nd, 



1 /2 

with Ld = {kBTe/q^ND) being the Debye screen- 
ing length associated with the semiconducting material. 
The parameter / depends only on the properties of the 
given semiconductor, while a depends also on the contact 
properties, thus being independent variables. Since the 
non-equilibrium current noise properties of devices with 
homogeneous stationary free carrier density profiles were 
shown to change qualitatively when passing from I < 1 
to I > 1, in what follows we will treat these two cases 
separately. 



(18) 



A. Homogeneous solution 

For q; = 1 [fic = No), and any value of I, there exists 
a trivial homogeneous solution to the model presented in 
Sec. Ill A, namely 



n{x) = Nd, E{x) = 



(19) 



Accordingly, the current- voltage (/ — V") characteristics 
is linear and satisfies Ohm's law / = V/Rbuik, with 
Rbuik = L/iqAfiNo) being the semiconductor resistance. 



B. Inhomogeneous solution 

For a > 1 (ric > Nd) the solution to the model pre- 
sented in Sec. Ill A is spatially inhomogeneous and can 
not be obtained in a closed analytical form. Only under 
asymptotic conditions it is possible to derive approximate 
analytical expressions^'^ (sec also appendix B). In any 
case, an exact solution of the transport model can be ob- 
tained by numerical integration following this procedure. 
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One first solves numerically Eq.(9) subject to the bound- 
ary conditions in Eq.(lO). Having found the electric field 
profile E{x), the carrier density profile n{x) is obtained 
from the Poisson equation (6). Furthermore, for a given 
value of the current /, the applied voltage is found as 
V = dxE{x), from where the I — V characteristic is 
obtained. 

In what follow we discuss the main features of the in- 
homogeneous solution. To this end we will consider two 
particular values of I, namely Z = 0.1 and Z = 50, as rep- 
resentative examples of the behavior observed for Z < 1 
and Z > 1, respectively. In both cases, the values of a 
will be varied in the range 1 < a < 10**, thus allowing us 
to explore both slightly and highly inhomogeneous situ- 
ations. 

Figure 1 displays the stationary free-carrier density 
profiles for I = 0.1 for several values of the applied bias 
and for a = 10 (Fig. la) and a = 10** (Fig. lb). As can 
be seen, the stationary free carrier density profiles are 
inhomogeneous and depend on the applied bias value. 

For a = 10 (Fig. la) the degree of inhomogeneity is 
rather small being at most of about 1%, and the free car- 
rier density departs more significantly from ric the lower 
the value of the applied voltage. The physical reason for 
this behavior can be found in the fact that the character- 
istic Debye screening length of the sample, roughly deter- 

1 /2 

mined by ric, i.e. Ldc = {ksTe/ q^n^ , is larger than 
the sample length L {1^ = L/Ldc = la^'^ = 0.1 x lO^/^ ^ 
0.3) thus not allowing the free carrier density to relax 
from its contact value Uc to its bulk value Nd, i.e. to the 
local charge neutral state. Similar behaviors appears for 
other values of Z < 1 and a > 1 as long as Ic = la^^'^ ^ 1. 

For a = 10^ (Fig.lb) the de gree of inhomogeneity be- 
comes appreciable being up to about a factor of 10. For 
low applied bias n{x) is almost symmetric and departs 
significantly from n^. The reason is that the character- 
istic Debye screening length, which can be roughly ap- 
proximated by Loc, is shorter than the sample length 
L {Ic = la'l'^ = 0.1 X 10^ = 10), thus allowing relax- 
ation towards the charge neutral state. However, due to 
the small value of l,.. complete relaxation of n(.x) to the 
value Ny) is not reached. At increasing applied bias, a net 
injection of carriers takes place resulting in an increase 
of the values of n{x) and in a rather asymmetric profile. 
Finally, for the highest applied biases carrier injection no 
longer takes place, and the density distribution is almost 
homogeneous with a value close to the contact value Uc- 
Similar behaviors are observed for other values of Z < 1 
and a > 1 as long as Ic = la^^"^ > 1. 

The current-voltage characteristics for I = 0.1 and for 
a ranging from 1 to 10^ are displayed in Fig. 2. For 
a < 10^, the I — V characteristics is found to remain lin- 
ear with a resistance given by Rc = L/{qA^nc)- This 
linear behavior can be understood by noting that for 
these values of a the free carrier density profile is quasi- 
homogeneous with a value approximately equal to Uc and 
almost independent of the applied bias (see Fig. la). 



Similar behaviors are found for other values of I < 1 and 
a > 1 as long as Ic = la^^^ < 1. For a = 10^ and iC the 
I — V characteristics in Fig. 2 displays a super-linear be- 
havior in an intermediate range of applied voltages, while 
it is linear at the lowest and highest applied bias. The 
resistance at the lowest bias, Riow , satisfies Rc < Riow < 
Rbuik, since in this range of bias the free carrier den- 
sity profile is almost independent of bias and satisfies 
ric < n{x) < No (see Fig.lb). Furthermore, the resis- 
tance for high applied bias Rhigh is almost equal to i?c, 
since in this range of bias the density profiles is almost 
independent of bias and given by n{x) ^ Uc (sec Fig. 
lb). Finally, the superlinear behavior is due to the net 
injection of carriers taking place in this bias regime (see 
Fig. lb) thus resulting in an increase of the conductivity 
of the sample. Similar behaviors appear for other values 
of Z < 1 and a > 1 as long as Ic = la^^^ > 1. 

Let us now consider the case of Z = 50. Figure 3 reports 
the stationary free carrier density profile for a = 10^ and 
several values of the applied bias. The qualitative behav- 
ior of the profiles is similar to those observed in Fig. lb, 
since both cases correspond to situations in which Zc > 1 
(in the present case Ic = la^^"^ = 50 x 10^/^ _ 1581). 
The main difference between these figures is that in Fig. 
3 relaxation to the local charge neutral state can take 
place for low bias, since the value of Ic is high enough. 
Similar behaviors are observed for other values of Z > 1 
and a > 1. 

The current- voltage characteristics for Z = 50 and sev- 
eral values of a are plotted in Fig. 4. The I — V charac- 
teristics is linear at the lowest bias, super-linear at inter- 
mediate bias and again linear at the highest bias. This 
behavior is similar to that reported in Fig. 2 for the 
curves with a = 10^ and 10"^, since both cases correspond 
to situations where Zc > 1. 

In the present case (Z > 1) one can reach an asymp- 
totic limit for Z ^ 1 and a ^ 1 in which the following 
analytical asymptotic expression for the I — V curve can 
be derived, (see appendix B) 



Hh 



1 qV 
I keT' 
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8 

I ksT^ 







- kBT ^ ' 



»J_ ( qV 

' \kBT 



/2 < 



qV 

ksT 

qV 

ksT 



(20) 



where wc have defined Iih = lkBT/{qRhuik)- Equation 
(20) gives a good description of the I — V curves for 
Z ^ 10 and a ^ 100, as can be seen in Fig. 4 where the 
symbols represent Eq. (20).^^ Note that in the asymp- 
totic limit, the low bias resistance is given by Riow = 
Rbuik, independent from any contact parameter since at 
low bias a state of local charge neutrality is achieved in 
the bulk of the sample (sec Fig. 3). Furthermore, in 
the intermediate bias regime, the I — V characteristics 
display a quadratic bias dependence, in agreement with 
the well-known Mott and Gurney law for space charge 
transport in diffusive conductors.'^ In this bias range, a 
strong injection of carriers from the contacts takes place, 
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as illustrated in Fig. 3. Finally, at the highest voltages 
the resistance is given by Rhigh = -Ro since the density 
profile is almost homogeneous and equal to ric (see Fig.3). 

V. CURRENT NOISE PROPERTIES 

We characterize the low frequency current fluctuation 
properties by means of the low frequency current spectral 
density defined as 

/+00 
{61{0)51{t)) dt. (21) 
-oo 

For the model presented in Sec. IIIB, Si{0) can be given 
a closed analytical expression that takes into account in 
an exact way the effects of the diffusion current on the 
non-equilibrium current fluctuations and hence is appli- 
cable to the whole range of system parameters. Below 
we discuss separately the cases corresponding to homo- 
geneous and inhomogeneous stationary density profiles, 
since they involve rather different mathematics. 



For / > 1 the noise properties change considerably. 
Accordingly, in the limit I ^ 1 Eq. (22) is well approxi- 
mated by^ 

r |forO<^<;i/3 

where S*j'^ = qlth- The above expression gives: Nyquist 
noise at low bias (0 < qV/ksT < l^^"^), a cubic 
dependence on current of noise at intermediate bias 
(Z'*/^ < qV/ksT < P) and shot noise at the highest 
bias {qV/ksT > Z^). The physical reason for the ap- 
pearance of shot noise in this case is the vanishing of 
long range Coulomb correlations because of a drift tran- 
sit time shorter than both the diffusion transit time and 
the dielectric relaxation time.'^ 



B. Inhomogeneous stationary profiles 



A. Homogeneous stationary profiles 

When the stationary profiles are homogeneous, Eq.(16) 
consists of a second order differential equation with con- 
stant coefficients. The solution of this equation, as well 
as an analytical expression for 5/(0), has been derived 
recently:'^ 

sm = ^ + 



Rbulk 

(Ai - A?) (e^i^ - 1) {e^-^ 
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1) (e^i^ - 1) - 



Ai (e^i^ + 1) (e^=^ 
where K = AqAkBTfrn, and 



1)] 



2Lr 



1 ± 




(23) 



with Le = ksTL/^qV). In Eq.(23) the subscript 1 corre- 
sponds to the plus term and the subscript 2 to the minus 
term. 

For ; < 1, Eq. (22) takes the form'^ 



Si{0) = 2qIcoth{qV/2kBT). 



(24) 



Equation (24) predicts a standard cross-over between 
Nyquist noise (Eq.(l)) and shot noise (Eq.(4)) for 
qV/ksT ~ 3. The physical reason for the appearance 
of shot noise is the absence of long range Coulomb corre- 
lations due to the fact that the sample length L is smaller 
than the Debye screening length Ld-^ 



When the stationary profile is inhomogeneous, Eq. 
(16) consists of a second-order stochastic differential 
equation with non-constant coefficients. To obtain an 
analytical solution of this equation, we use a method de- 
veloped in Refs. [ 13-15]. The method is based on the 
fact that 



p(x) = dE{x)/dx, 

constitutes a particular solution of Eq. (16). 
basis one can find an analytical expression for 
trie field fluctuations 5E^.{t) that solves Eq. 
satisfies the boundary condition in Eq. (17).^^ 
expression of the electric field fluctuations we 
the voltage fluctuations, SV{t) = J^^ SEx{t)dx. 
ing a procedure similar to that outlined in Ref. 
voltage fluctuations can be expressed as: 



with 



5V{t) 



eAD 



f 
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WZ{x)6I^{t)dx, 



(26) 
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Here, (f>{x) is the electrie potential. The function VZ{x) 
is referred to as the impedance field, since it satis- 
fies 



Z{0) 



V Z(x)dx 



.-F X -1 



dJ_ 
dV 



(30) 



From Eq. (27), the low frequency spectral density of the 
voltage fluctuations is given by 

S'y(O) = AAq^D [ [VZ{x)]^n{x)dx, (31) 
Jo 

where use is made of Eq. (12) and (13). The low fre- 
quency spectral density of current fluctuations is then 
evaluated from^° 



5,(0) 



MO) 

m 



2 ' 



(32) 



Equations (28)- (32) constitute an exact solution of the 
noise model presented in Sec. Ill for the general case 
of an inhomogeneous stationary profile. They represent 
the main result of this paper. Note, that the solution 
incorporates the effects of the diffusion current without 
any approximation, thus allowing us to investigate the 
transition from homogeneous to highly inhomogeneous 
stationary profiles. 

In order to evaluate Eqs. (28)-(32) we will use the 
steady-state profiles calculated in Sec. IV B. As for the 
case of transport we will consider the values I = 0.1 and 
I = 50, as representative examples of the behavior ob- 
served for / < 1 and ^ > 1. 

Figure 5 reports the low frequency spectral density of 
current fluctuations 5/(0) as a function of the current 
as obtained from Eqs. (28)-(32) for I = 0.1 and several 
values of a. For comparison the results for the homo- 
geneous case, as calculated from Eq. (22), are also dis- 
played (dashed line). We note that for the lowest degree 
of inhomogeneity (a < 10^) the non-equilibrium current 
noise results can be well approximated by Eq. (24) (sym- 
bols in Fig. 5). The reason is that for these values of a 
the profiles are quasi-homogeneous with free carrier den- 
sity Uc and with Debyc screening length satisfying Ic < 1, 
which are the conditions for the validity of Eq.(24). Sim- 
ilar behaviors are observed for other values of Z < 1 and 
a > 1 for which Ic < 1. 

By contrast, for a high degree of inhomogeneity {a — 
10'^ and 10^) the results deviate from Eq.(24), and re- 
gions where the current noise is suppressed below the 
shot noise limit are observed at intermediate bias values. 
The bias values where shot noise siipprcssion takes place 
are found to correspond to those applied bias for which 
significant net carrier injection from the contacts takes 
place (see Fig. lb) . The physical origin of the shot noise 
suppression should be traced back to the correlations in- 
duced by the long range Coulomb interaction which are 
active in this range of bias. Similar behaviors are ob- 
served for other values of / < 1 and a > 1 for which 
Ic > 1. 



Figure 6 reports the low frequency spectral density of 
current fluctuations 5/(0) as a function of the current as 
obtained from Eqs. (28)-(32) when I = 50 and for several 
values of a. For comparison, the results for the homoge- 
neous case, as calculated from Eq.(22), are also displayed 
(dashed line). As can be seen, a qualitative new behavior 
at intermediate current values is identified with respect 
to the homogeneous behavior for all values of a. The new 
behavior is found to tend to a square root dependence on 
the current and to cover broader current intervals as a is 
increased. This fact can be quantitatively evaluated by 
means of the following explicit expression valid for Z ^ 10 
and a > 100 
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(see Fig. 6 where the symbols correspond to the ap- 
proximate expression). Otherwise, for 1 < Z < 10 and 
1 < a < 100 the exact result should be used. 

The dependence of the noise power as the square root 
of the current corresponds to the well known double ther- 
mal noise regime found in space charge limited deviccs^'^ 
(see Eq.(2)). We note that the double thermal noise 
behavior is restricted to current values in the range 

Z/18 < I /Ith ^ (4.5q:*Z)"^^''. For higher currents the noise 
properties deviate from the double thermal noise behav- 
ior and start resembling those corresponding to homoge- 
neous conditions (cubic current dependence followed by 
shot noise). This fact is illustrated in Fig. 7 where we 
compare the exact results of the present paper (continu- 
ous line) with those of existing theories for double ther- 
mal noise represented by Eq. (3) (triangles) for I = 50 
and a = 10"^. As can be seen, by neglecting the diffusion 
current, existing theories can only be applied up to bias 
values below the onset of the cubic region. At applied 
bias above this onset, the diffusion current plays a rel- 
evant role and the exact theory presented here must be 
used. It is worth remarking that for devices operating 
under strong current injection conditions (a > 10"* ) the 
quasi-homogenous behavior predicted at the highest bias 
can be hardly observed in practice because of hot-carrier 
effects. However, in the low current injection regime 
(1 < a < 10"^) this behavior should be experimentally 
accessible as discussed in the following suggested exam- 
ple. 



C. Example 



As an example to illustrate the experimental acces- 
sibility of the theoretical predictions described above 
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wc consider the following particular case. As semicon- 
ductor material we consider high resistive p-type CdTc 
at room temperature with a free carrier density (holes) 
p = 10^ cm~'^. These low values of the free carrier den- 
sity can be obtained by means of compensation.^^ The 
hole mobility is taken to be /z ~ 40 cm^/Vs, the effec- 
tive hole mass m* =0.8 mo, where toq is the electron 
mass, and the dielectric constant e = 10.3 eo, with eo 
being the vacuum permittivity. The sample length is 
assumed to be L = 4 mm and the cross sectional area 
^ = 40 mm^. As metal for the contacts we consider 
gold, which has been reported to form almost ideal metal- 
semiconductor junctions on CdTe.^^ The value for the 
Au/CdTe barrier height for holes, deduced from that for 
electrons, is = 0.55 V, where use has been made 
of the relationship 4'bp = Eg — (ptn with Eg = 1.48 eV 
and (j>bn = 0.93 V.^^ According to the value of (pbp the 
contact density for holes is Pc = 10^" cm~^, where we 
used that Pc = Nvexp[—q(l)bp/{kBT)], with Ny = 1.8 
10"'^'' cm^'"' being the CdTe effective density of the states 
in the valence band at room temperature. The metal 
contact described above forms an ohmic injecting con- 
tact in the terms described in the present paper (see also 
at the end of appendix A.) 

For the set of parameters considered above, we have 
Lo = 0.4 mm and Ldc = 0.04 mm, from where I = 
L/Ld = 10, a = Uc/Nd = 100 and 1^ = L/Ldc = 100. 
Figure 8 displays the calculated current / (right axis) 
and low frequency current spectral density (left axis) as 
a function of the applied. According to the calculations, 
the predicted I — V characteristic is linear up to bias 
voltages around 0.6 V. For higher bias it tends to be 
quadratic, up to F ~ 300 V, from where it returns to lin- 
earity. Furthermore, the calculated current spectral den- 
sity S'/(0) displays the Nyquist thermal value up to volt- 
ages around 0.3 V. Then, it increases with voltage accord- 
ing to the double thermal noise behavior up to around 20 
V. At further increasing voltages, 5*7(0) increases sharply 
with voltage according to the cubic dependence with ap- 
plied current up to y ~ \kV , where shot noise appears. 
It is worth noting that the calculated values of the low fre- 
quency current spectral density are well inside the range 
of experimental accessibility (state of the art correlation 
spectrum analyzers^" are able to reach noise levels as 
low as 10~^^ A^/Hz). Moreover, the electric fields reach 
maximum values up to Eav = V/L = 5 kV/cm, which 
are still below that for the onset of hot electron effects in 
p-type CdTe. We conclude that such an example shows 
that the theoretical predictions presented in this paper 
are accessible to an experimental confirmation. 



VI. CONCLUSIONS 

Wc have carried out an analytical theory of trans- 
port and current fluctuation properties in metal- 
semiconductor-metal structures made of highly resistive 



semiconductors. The theory includes the effects of the 
diffusion current in an exact way, thus allowing us to 
study the whole range of physical conditions concern- 
ing the strength of the applied bias and the level of car- 
rier injection from the contacts. It is shown that in the 
low carrier injection limit for which Ic = L/Ldc < 1, 
where L is the sample length and Ljjc the Debye screen- 
ing length associated with the free carrier density at the 
metal-semiconductor interface ric, the structure behaves 
like a linear resistor with low frequency noise proper- 
ties given by 5/(0) = 2ql coth[qV/ {2k bT)]. In the inter- 
mediate carrier injection regime, roughly determined by 
1 < ^ 100, the structure displays a linear-superlinear- 
linear current-voltage characteristics. In this regime, 
the current spectral density displays a cross-over from 
Nyquist noise to shot noise mediated by a region depend- 
ing first as the square root of the current (double ther- 
mal noise) and then as the third power of the current. 
Finally, under strong carrier injection conditions Zc ^ 1 
the standard theory of space charge limited diodes is re- 
covered. In this limit the current-voltage characteristics 
is first linear and then quadratic. Accordingly, the cur- 
rent spectral density displays Nyquist thermal noise at 
low bias followed by double thermal noise at higher ap- 
plied bias. We suggest that high resistive p-type CdTe 
is one of the best suited materials to provide an experi- 
mental test of the theoretical predictions in the small and 
moderated injection regimes. 
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APPENDIX A: METAL-SEMICONDUCTOR 
MODEL CONTACT 

In this appendix we justify the boundary conditions 
used along the paper for both the transport, Eq. (7) and 
the noise properties, Eq. (15). Following Refs. [ 1, 22, 
23], the general boundary conditions for the transport 
through a metal-semiconductor contact read 



n(0) 



I+Jl 
qAvO ' 



n{L) 



-I + Ij 
qAv^ 



(Al) 



with and (res. and I^) being the recombination 

velocity and saturation current, respectively, of the con- 
tact located at x = (res. x = L). The currents 1^ and 
are given by = qAv^n^ and 1^ = qAv^rf^, where 
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n, 



eg 



,69 



Nce~ 



(A2) 



Here, Nq is the effective density of states in the con- 
duction band and (j)^^ (res. is the barrier height at 
contact a; = (res. x = L). For the case of ideal metal- 
semiconductor junctions, the barrier heights are given by 
'Pbn = 4>°m- X and (/)^„ = 4,^- X, with x being the 
semiconductor affinity and (f)^ (res. 0^) the work func- 
tion of the metal located at x = (res. x = L). For 
highly resistive materials the diffusion approximation^ 
can be applied. In this approximation it is assumed 
that the current takes values much smaller that the 
thermionic current, i.e. / <C This condition implies^^ 
T/hh < Ts/hh = O.II3, where 



(3 = ^i 



2TTNDm* 



a 



Nd' 



(A3) 



In this limit, one can approximate Eq.(Al) by 



x/Ld, L L/Ld, U{x) U{x)/kBT, where Eth = 

1 /2 

(fcBTiV£,/e)'/^Equation (Bl) can be easily integrate be- 
tween a; = and x = L to give the following relation 



EL + nn[l-^ \ -EQ + I\n[l 



^]=L. (B2) 



where El = E{L) and Eq = E{0). On the other hand, 
by using again Eq.(Bl), one can derive an explicit ex- 
pression for the potential energy difference U{x) — Ul = 
E{x)dx, which once evaluated at a; = gives for the 
applied bias, V = Ul — Uq, 



2 

V =^+IEL+tln 



-2 



El 
J 



(B3) 



n(0) 



qAv° 



eq ^ 
^0 ' 



n{L) 



qAv^ 



(A4) 



which for the case of a symmetric structure correspond 
to the boundary conditions in Eq.(7). This result also 
implies that the resistance of the metal-semiconductor 
interface is negligible small in comparison to the bulk 
resistance. As a result, the fluctuations generated at 
the metal-semiconductor interface can be completely ne- 
glected. This is equivalent to approximate^^ 



Sno{t) ~ 0; SnL{t) ^ 0, 



(A5) 



in agreement with the boundary conditions used along 
the paper Eq. (15). 

For the concrete example considered in Sec.VC, one 
has P = 9 10~^ ^ 1, which allows us to apply the dif- 
fusive approximation up to currents satisfying the con- 
dition Ig/Ith = ct/P = 1-1 10''', well above the current 
values required to observe the different behaviors here 
identified. 



To end up with an explicit expression for the I — V char- 
acteristics, we impose the boundary condition El = I/a 
in Eqs.(B3) and (B2), thus giving the I — V characteris- 
tics in parametric form 



V 



L2 [i+in(l-i) -no-ln(l-«o) 
L 

[i + ln(l-i)-t.o-ln(l-«o)]' 



(B4) 



(B5) 



where the parameter uq = Eq/I satisfies the condition 
1/a < Uq < 1. By further expanding the previous ex- 
pression for a ;» 1 one then arrives at Eq. (20). 



APPENDIX B: ASYMPTOTIC TRANSPORT 
THEORY 

In this appendix we derive the asymptotic expression 
for the I — V characteristics used in Eq.(20) of Sec. IV B. 
To this purpose, we start by considering the limit when 
1^1. In this limit, one can neglect the diffusive contri- 
bution in Eq.(9), thus obtaining a first order differential 
equation for the electric field of the form 



-E{x) 



dE{x) 
dx 



1 



(Bl) 



subject to the boundary condition E{L) = I /a. In 
equation (Bl), and along this appendix, we use di- 
mensionless variables to simplify the expressions. Ac- 
cordingly, we take E —>■ E/Eth, I — > I/Ith, x —>■ 
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FIG. 1. Free carrier density profiles normalized to the 
donor density for several values of the applied bias, for Z = 0.1 
and (a) a = 10 (Zc = 0.3) and (b) a = lO" {Ic = 10). 
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FIG. 2. Current voltage characteristics for I — 0.1 and sev- 
eral values of a (the corresponding value of Ic is displayed for 
each curve). The dashed line corresponds to the homogeneous 
solution. 
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FIG. 4. Current voltage characteristics for Z = 50 and sev- 
eral values of a (the corresponding value of Ic is displayed 
for each curve). The dashed line corresponds to the homo- 
geneous solution. The symbols correspond to the asymptotic 
expression of the I — V characteristics given in Eq. (20). 
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FIG. 5. Low frequency current spectral density normalized 
to Sith. = Qlth as a function of the current normalized to Ith 
for I — 0.1 and several values of a (the corresponding value of 
Ic is displayed for each curve). The dashed line corresponds 
to the homogeneous solution. The symbols correspond to the 
coth-like expression in Eq. (24). 
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FIG. 6. Low frequency current spectral density normalized 
to Sith = ihh as a function of the current normalized to Ith 
for I = 50 and several values of a (the corresponding value of 
Ic is displayed for each curve). The dashed line corresponds 
to the homogeneous solution. The symbols correspond to the 
approximate expression in Eq. (33). 
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FIG. 7. Electric current (dashed line) and low frequency 
current spectral density (solid line) as a function of applied 
bias for I = 50 and a = 10^. The different behaviours are 
schematically indicated. The current is normalized to Ith and 
the spectral density to Sj^ = qltn- For comparison the results 
obtained from Eq. (3) are also displayed as triangles. 
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FIG. 8. Electric current (dashed line) and low fre- 
quency current spectral density (solid line) as a func- 
tion of applied bias for a highly resistive p— type CdTe 
metal-semiconductor-metal structure. Parameters: free car- 
rier density 10* cm"'', sample length 4 mm, cross sectional 
area 40 mm^, temperature 300 K and contact density 10^^° 
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